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Pluripotent cells possess the ability to differentiate into 
any cell type. Commitment to differentiate into specific 
lineages requires strict control of gene expression to 
coordinate the downregulation of lineage inappropriate 
genes while enabling the expression of lineage-specific 
genes. The nucleosome remodelling and deacetylation 
complex (NuRD) is required for lineage commitment of 
pluripotent cells; however, the mechanism through which 
it exerts this effect has not been defined. Here, we show 
that histone deacetylation by NuRD specifies recruitment 
for Polycomb Repressive Complex 2 (PRC2) in embryonic 
stem (ES) cells. NuRD-mediated deacetylation of histone 
H3K27 enables PRC2 recruitment and subsequent H3K27 
trimethylation at NuRD target promoters. We propose a 
gene-specific mechanism for modulating expression of 
transcriptionally poised genes whereby NuRD controls 
the balance between acetylation and methylation of his- 
tones, thereby precisely directing the expression of genes 
critical for embryonic development. 
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Introduction 

Precise control of gene expression is essential both for 
embryonic stem (ES) cell differentiation and for early em- 
bryonic development. Lineage-specific genes must be 
expressed in a precise temporal and spatial fashion, but it is 
equally important that expression of other genes is repressed. 
In short, it is the balance between transcription of lineage 
appropriate genes and repression of inappropriate genes that 
allows progression through development. 

Regulation of gene expression is intimately linked to 
chromatin state, which in turn is heavily influenced by the 
presence of post-translational modifications of the histone 
proteins contained within the nucleosome. These modifica- 
tions include acetylation, methylation, phosphorylation, 
ubiquitylation and sumoylation at specific residues both on 
the histone globular domain and on the histone N-terminal 
tails (Berger, 2007; Kouzarides, 2007). The relationship 
between specific histone modifications and transcriptional 
state has been well established. For example, acetylation of 
histone tails is generally associated with active transcription, 
while methylation may be linked to either activation or 
silencing of transcription depending on which histone residue 
is modified (Jenuwein and Allis, 2001; Rice and Allis, 2001; 
Roh et al, 2005; Kouzarides, 2007; Shahbazian and Grunstein, 
2007). In some cases, such as at H3K27 and H3K9, acetylation 
or methylation can occur at the same histone residue and it is 
the balance between these opposing modifications that 
determines the transcriptional status of that region (Rice 
and Allis, 2001; Tie et al, 2009; Jung et al, 2010; Pasini 
et al, 2010b). 

The nucleosome remodelling and deacetylation (NuRD) 
complex is a transcriptional co-repressor essential for devel- 
opmental transitions in early embryogenesis as well as for ES 
cell function (reviewed in McDonel et al, 2009). In the 
absence oiMbd.3, which encodes a core structural component 
of NuRD, the complex does not form and embryonic devel- 
opment stalls at the implantation stage, with the mutant 
embryos failing to form differentiated cell types (Zhang 
et al, 1999; Hendrich et al, 2001; Kaji et al, 2006, 2007). ES 
cells lacking NuRD are viable but are unable to exit self- 
renewal and commit to differentiation upon withdrawal of 
LIF (Kaji et al, 2006) . NuRD components have been reported 
to interact with Oct4, a protein essential for the maintenance 
of pluripotency in ES cells (Liang et al, 2008; Pardo et al, 
2010; van den Berg et al, 2010), despite the fact that several 
NuRD components are dispensable for pluripotency 
(McDonel et al, 2009). NuRD function is also important for 
homeostasis of both haematopoietic and epithelial stem cells 
(Williams et al, 2004; Kashiwagi et al, 2007) . Aberrant gene 
expression patterns have been demonstrated in embryonic 
and somatic cell types in the absence of a functional NuRD 
complex (Kaji et al, 2006, 2007; Yoshida et al, 2008). Taken 
together, these findings suggest that NuRD-mediated gene 



©2012 European Molecular Biology Organization 



The EMBO Journal VOL 31 | NO 3 | 2012 593 



NuRD and PRC2 interact to control gene expression 

N Reynolds et al 



regulation is required for stem cell fate decisions both in 
culture and during embryonic development. 

The NuRD complex has been extensively characterised 
biochemically and contains multiple protein subunits, includ- 
ing the class I histone deacetylases HdacI and II, and the ATP- 
dependent chromatin remodelling component Mi2|3 (re- 
viewed in McDonel et al, 2009). Nevertheless, the precise 
mechanism through which NuRD controls gene expression is 
unclear. While a reduction of histone H3K27 trimethylation 
was found in plants lacking the NuRD component PICKLE 
(Zhang et al, 2008a) and knockdown of Mi2 in mammalian 
cells was shown to be associated with a decrease of H3K27 
methylation at one target locus (Morey et al, 2008), no global 
changes in histone modifications have been detected upon 
loss of NuRD in mammalian cells (Kaji et al, 2006). As yet, no 
direct link between the action of this histone deacetylase 
containing complex and the methylation state of H3K27 has 
been demonstrated. 

Like NuRD, Polycomb group (PcG) complex proteins re- 
press transcription in ES cells and during multiple develop- 
mental programs. Biochemically, Polycomb-Repressive 
Complex 2 (PRC2) acts to di- and tri-methylate H3K27 via 
the methyltransferase activity of Ezh2 (Cao et al, 2002; 
Czermin et al, 2002; Kuzmichev et al, 2002; Muller et al, 
2002). This may then lead to the recruitment of PRC1 which 
ubiquitinates histone H2A, thereby silencing transcription of 
target genes (Muller and Verrijzer, 2009; Christophersen and 
Helin, 2010). In the absence of PcG proteins, both early 
embryonic development and ES cell differentiation are dis- 
rupted, although ES cells remain pluripotent (Faust et al, 
1995; O'Carroll et al, 2001; Voncken et al, 2003; Pasini et al, 
2004; Isono et al, 2005; Boyer et al, 2006). Polycomb targets 
have been extensively identified and include genes with 
'bivalent' modifications, that is, a combination of H3K4 
trimethylation, associated with active transcription, and 
H3K27 trimethylation, a repressive transcriptional mark. It has 
been proposed that bivalent genes in ES cells are poised for 
transcription and may be activated or silenced by the removal of 
the repressive or active marks, respectively (Azuara et al, 2006; 
Bernstein et al, 2006; Mikkelsen et al, 2007). 

While there is some evidence for an in-vivo interaction 
between PcG proteins and the NuRD complex in various 
organisms (Kehle et al, 1998; Unhavaithaya et al, 2002; 
Morey et al, 2008; Aichinger et al, 2009), the precise nature 
of this interaction has not been characterised. Nevertheless, 
the importance of a balance between the acetylation and 
methylation state of H3K27 has been shown in both 
mammalian cells and flies (Tie et al, 2009; Jung et al, 2010; 
Pasini et al, 2010b) and could provide the link between these 
two complexes in stem cell function. 

By comparing levels of specific chromatin modifications in 
ES cells with or without functional NuRD complex, we 
demonstrate the role played by NuRD in regulating the 
balance between acetylation and methylation state of 
H3K27. We propose a two-step model for repression of gene 
expression through the combined action of the deacetylase 
activity of NuRD and the methlytransferase activity of PRC2, 
which provides a molecular mechanism underlying 
NuRD-mediated lineage commitment of ES cells. Moreover, 
we illustrate a new level of complexity in transcriptional 
regulation through PcG proteins by showing that PRC2 can 
be directed to act at specific genes by NuRD. 

594 The EMBO Journal VOL 31 | NO 3 | 2012 



Results 

Gene expression changes in absence of Mbd3/NuRD in 
ES cells 

To obtain a global view of NuRD-mediated transcriptional 
regulation, gene expression profiles of wild-type and Mbd3- 
null ES cells were compared by microarray analysis. 
A surprisingly large number of genes was found to be 
significantly downregulated (839 genes, P<0.01) while 531 
genes showed significant derepression (P<0.01) in Mbd3~ /_ 
compared with wild-type ES cells (Supplementary Table 1). 
While many of these changes in gene expression will be 
caused indirectly by a lack of NuRD, these numbers indicate 
that NuRD is likely to activate as well as silence gene 
transcription in ES cells. To investigate the means by which 
NuRD acts to repress transcription, we decided to focus on 
the upregulated genes identified in this study. Expression 
changes were verified for a subset of upregulated and control 
genes by quantitative RT-PCR in both Mbd.3^ and Mbd3 fIOX/ ~ 
ES cells (Figure 1A). 

Although NuRD has long been known to be a transcrip- 
tional silencer, only 0.2% of the genes upregulated in 
Mbd3^^ ES cells show hallmarks of transcriptionally inac- 
tive chromatin, such as H3K27 trimethylation, in wild-type 
ES cells (according to Mikkelsen et al, 2007). In contrast, 17% 
of upregulated genes in wild-type cells are associated with 
'bivalent' chromatin, while 64% are associated with 
H3K4me3 but not with H3K27me3. The proportions of 
genes with either H3K4me3 alone, or with both H3K4me3 
and H3K27me3 that are upregulated in the absence of NuRD 
are similar to those seen on a genome-wide scale in ES cells 
(Mikkelsen et al, 2007). This indicates a lack of specificity 
towards these particular histone modifications at loci targeted 
by NuRD. While the bivalent mark has been associated with 
poised genes, H3K4me3 marks active genes (Kouzarides, 
2007; Stock et al, 2007). Therefore, rather than functioning 
as a transcriptional silencer, our analysis indicates that NuRD 
acts to modulate the output of both 'poised' and actively 
transcribed genes in ES cells. 



NuRD binds directly to loci with H3K4me3 and 
H3K4me3/H3K27me3 

Chromatin immunoprecipitation (ChIP) was carried out using 
an antibody specific to Mi2p\ a defining component of the 
NuRD complex, to determine which of the misregulated 
genes were direct targets of NuRD (Figure IB). ChIP profiles 
for Mi2f3 at physiologically relevant target genes using this 
antibody are highly similar as those for a tagged protein 
(unpublished observations). In addition, ChIP using antibo- 
dies to other components of the NuRD complex, Mta2 and 
Hdacl, as well as an Avi-tagged Mbd3 showed comparable 
binding properties to Mi2p (Supplementary Figure SI). We 
have seen that Mi2p is generally associated with chromatin, 
but that regions of specific enrichment can be identified (NR 
and BH, unpublished observations). Therefore, we focussed 
on regions of the genome close to the transcription start sites 
of differentially expressed genes. Genes at which Mi2|3 was 
found to be enriched include those with both H3K27me3/ 
H3K4me3 (e.g., Ppp2r2c, Html, Klf4 and Tbx3) and 
H3K4me3 only-associated promoters (e.g., Smad7, Sohlh2, 
Klf2, KlfS and Lefty 2) (Figure IB). 
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Figure 1 Identification of direct gene targets for NuRD and PRC2. [A) Quantitative RT-PCR comparing transcript levels inMfod3~ / ~ ES cells to 
those in wild-type ES cells. Results are plotted as log 10 fold change relative to wild-type levels. Error bars indicate standard error of the mean 
(s.e.m.). (B) Chromatin IP in wild-type cells for either Mi2[3 or IgG control, with qPCR for proximal promoter regions of genes shown. Results 
are plotted as percentage of input DNA. Histone signatures (according to Mikkelsen et al, 2007) are indicated underneath. Asterisks denote loci 
at which ChIP for Mi2p is significant with respect to IgG control fP< 0.005) . (C) Chromatin IP in wild-type cells for either Suzl2 or IgG control, 
plotted as percentage of input DNA. Histone signatures are indicated underneath. Asterisks denote loci where ChIP for Suzl2 is significant with 
respect to IgG control (P< 0.005). Error bars indicate s.e.m. 



To obtain a more global view of NuRD-associated regions 
in ES cells, we performed ChIP followed by high-throughput 
sequencing (ChlP-seq) for Mi2p in wild-type ES cells. The 
results obtained are in good agreement with our ChlP-qPCR 
data; however, the affinity of antibodies to native Mi2|3 is 
relatively low and the high-throughput data are inexhaustive 
(Supplementary Table 2) . 
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PRC2, like NuRD, is a transcriptional repressor complex 
required for gene silencing in ES cells. PRC2 acts by directing 
methyltransferase activity to di- and tri-methylate H3K27. 
ChIP assays performed with an antibody to Suzl2, a core 
component of the PRC2 complex, confirmed that binding of 
PRC2 is restricted to regions that are enriched in H3K27me3 
(Figure 1C). NuRD components were found at a subset of 
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Figure 2 Comparison of multiple histone modifications between Mb63~'~ and wild-type ES cells. (A) ChIP for histone modifications. 
Enrichment at promoter regions of genes for histone modifications, bulk histone H3 levels or IgG control plotted as null value relative to wild 
type. Error bars indicate s.e.m. (B) Whole histone extracts from wild-type, Mbd3 _/ ~ and Eed mutant ES cells blotted with antibodies for histone 
H3 (loading control), H3K27me3 and H3K27ac to show relative levels. 



bivalent promoters as well as those with only H3K4 trimethy- 
lation, revealing that while some overlap is observed in target 
loci between NuRD and PRC2, the two complexes target 
distinct sets of genes (Figure IB and C). 

H3K27 is the focus of histone modification by NuRD 

NuRD is a chromatin remodelling complex exhibiting 
deacetylase activity. We therefore, expected that loss of reg- 
ulation of gene expression in the absence of NuRD would be 
reflected in changes to chromatin state, particularly in histone 
acetylation levels. To test this hypothesis, we used ChIP to 
compare levels of various histone modifications between 
Mbd3 flox/ ~ and Mb<33~'~ cells for four NuRD target genes 
and two non-targets, all associated with bivalent chromatin 
domains. Promoter-proximal regions were assayed by qPCR 
from chromatin purified with antibodies specific to either 
acetylation or methylation of histones H3 and H4, and fold 
change determined relative to wild type (Figure 2A). In 
general, there was a strong association between an increase 
in levels of marks of active transcription (e.g., histone acet- 
ylation), a corresponding decrease in repressive marks 
(H3K27me3 and H3K9me3) and increased transcription levels 
in MM3-null cells (Figures 1A and 2A). In contrast, no clear 
pattern in changes to H3K4me3 levels was apparent between 
Mbd3 flox/ ~ and Mbd3 _/ ~ ES cells. Importantly, comparison of 
ChIP levels for the same regions using an antibody specific to 
histone H3 showed no overall change, indicating that the 
differences detected for specific histone modifications were 
not due to variation in nucleosome occupancy. 

While a correlation between an increase in transcription 
and an increase in chromatin marks of active transcription 
was not unexpected, the most profound changes seen in 
Mbd3~ / ~ ES cells are the increase in H3K27 acetylation and 
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loss of H3K27 trimethylation (Figure 2A). These changes are 
localised specifically to NuRD target genes, but not to the 
non-targets assayed (Figure 2A) . A similar effect was seen for 
H3K9 acetylation and trimethylation, although the magnitude 
of changes observed was less striking. These reciprocal 
changes in methylation and acetylation at H3K27, and to a 
lesser extent at H3K9, in the absence of Mbd3 led us to 
hypothesise that NuRD acts through deacetylation of these 
residues to control chromatin state and thereby the transcrip- 
tional status of specific genes. 

PRC2 is required to maintain global H3K27 methylation 
levels in ES cells (Montgomery et al, 2005). Similarly, we find 
that NuRD is required to maintain H3K27Me3 levels at 
targeted bivalent genes (Figure 2A). To determine whether 
NuRD is also required for global H3K27 trimethylation levels, 
we next compared H3K27 acetylation and trimethylation 
levels in bulk histones prepared from Mbd3T , ~ and £ed _/ ~ 
ES cells. In the absence of Eed, we observed a genome-wide 
loss of H3K27 trimethylation and an increase in acetylation 
levels (Figure 2B), consistent with published reports (Tie 
et al, 2009; Pasini et al, 2010b). In contrast, no changes in 
global levels of either modification were detectable by 
western blot in bulk histones extracted from Mbd3-null 
lines, consistent with our hypothesis that NuRD controls 
H3K27 modifications only at specific target genes. 

Reciprocal changes in H3K27 acetylation and 
methylation are common to bivalent NuRD target loci 

To address the extent to which this switch between acetyla- 
tion and methylation status at H3K27 correlates with the 
action of NuRD in ES cells, we performed ChlP-seq for 
H3K27ac and H3K27me3 in parental and Mbd3~'~ ES cells. 
Full details of genomic regions associated with each modifi- 

©2012 European Molecular Biology Organization 



NuRD and PRC2 interact to control gene expression 

N Reynolds et al 



cation in wild-type and null cells are listed in Supplementary 
Table 3. 

Comparing average ChlP-seq profiles for wild-type sam- 
ples, bivalent genes displaying altered expression in the 
absence of Mbd3 also showed a propensity for increased 
trimethylation and decreased acetylation at H3K27 relative 
to the average levels for all genes in the sample (Figure 3A; 
Supplementary Figure S2; Supplementary Table 3). In the 
absence of NuRD, differentially expressed bivalent genes 
showed a loss of H3K27 trimethylation, particularly around 



the transcription start site, combined with a reciprocal in- 
crease in the relative abundance of H3K27 acetylation 
(Figure 3A). While these changes in H3K27 acetylation and 
trimethylation are apparent across the entire gene body and 
at the transcription termination site (Supplementary Figure 
S2), they are most pronounced close to the transcription start 
site (Figure 3A). The same effect was seen when the compar- 
ison was repeated using Mi2p targets as determined by ChlP- 
seq (Supplementary Table 2; Supplementary Figure S3). Here, 
an overall decrease in trimethylation coincident with an 
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increase in acetylation at H3K27 was evident. These results 
are in agreement with the hypothesis that these two histone 
marks act in opposition and that their relative levels are 
controlled by NuRD activity. 

We sought to confirm the patterns revealed by the ChlP- 
seq data at specific, bivalent loci for a NuRD target {Html) 
and a non-target (T, also known as Bmchyury) . T expression 
was not altered in Mbd3-null cells, whereas Html increases 
by ~ 10-fold in null cells (Figure 1A; Supplementary Table 1). 
Detailed binding profiles for aH3K27ac, cdT3K27me3 and 
ocMi2p were compared across a ~6Kb region spanning the 
transcription start sites for these genes in both Mbd3 flox/ - and 
Mbd3T'~ ES cells (Figure 3B). Mi2p ChIP supports the 
hypothesis that NuRD binds across the Html promoter region 
and within the body of the gene, but is specifically enriched 
close to the transcription start site. Mi2|3, however, is present 
at very low levels across the equivalent region of the T 
promoter (Figure 3BJ. H3K27 trimethylation is detected 
across both the Html and T loci in wild-type cells, while 
very little H3K27 acetylation is present at either locus. In 
contrast, ES cells lacking Mbd3 display an absence of H3K27 
trimethylation and a significant increase in H3K27 acetylation 
at Html. These changes are not evident at the T locus, 
suggesting that the effect is specific to NuRD target genes. 
At Html, a strong peak of Mi2p binding is observed close to 
the transcription start site that coincides with the peak of 
H3K27 acetylation present in Mbd3 _/ ~ cells. Thus, we 
conclude that NuRD is targeted to a subset of bivalent 
genes where it is required to prevent H3K27 acetylation 
and, indirectly, to maintain H3K27 methylation. 



Gene-specific recruitment of PRC2 is dependent on 
NuRD 

Thus far, our data indicate that NuRD and PRC2 work in 
concert at bivalent NuRD target promoters to maintain H3K27 
in a deacetylated, trimethylated state. We next addressed the 
nature of the interplay between NuRD and PRC2 at their target 
genes. We were unable to detect a direct protein-protein 
interaction between the two complexes in ES cells by immu- 
noprecipitation (Supplementary Figure S4), suggesting that 
neither complex is physically recruited by the other to specific 
chromatin sites. Further, levels of PRC2 component proteins 
were essentially unchanged in Mbd3~ / ~ ES cells, and NuRD 
subunits were present at approximately wild-type levels in 
Eed~'~ ES cells (Figure 4A). These results do not support an 
interdependence at the level of component protein abundance 
or stability between the two complexes (Figure 4A) . 

We next asked whether NuRD activity was required for 
PRC2 to be physically recruited to its target loci. ChIP for 
Suzl2 in Mbd3 flox/ ~ and Mbd3^ cells revealed a loss of 
PRC2 at NuRD target genes in the absence of NuRD, whereas 
Mi2p levels were not significantly changed at its target loci in 
Eed mutant ES cells (Figure 4B and C) . ChIP for Jarid2, which 
directs PRC2 to target loci (Peng et al, 2009; Shen et al, 2009; 
Landeira et al, 2010; Pasini et al, 2010a) also revealed a 
dependency on NuRD for its recruitment while no change 
was seen in Mi2p binding in Jarid2-null ES cells (Shen et al, 
2009; Figure 4D and E). This observation supports a model in 
which recruitment of PRC2 at NuRD target loci is dependent 
on either the presence or activity of NuRD, but in which the 
loss of PRC2 has no effect on Mi2p recruitment. 
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Figure 4 Interdependency between NuRD and PRC2 chromatin interactions. (A) Western blot showing relative levels of NuRD and PRC2 
components in wild-type, Mbd3 _/_ and Eed mutant cells. Alpha tubulin is used as loading control. (B) ChIP for Suzl2 in Mfad3-null cells 
relative to wild type. (C) ChIP for Mi2p in Eed mutant cells relative to wild type. Association of Mi2p and Suzl2 in wild-type cells at each gene is 
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type. ChIP data shown is for bivalent genes only. Error bars indicate s.e.m. 
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We next sought to determine whether it is the physical 
presence or deacetylase activity of NuRD that is responsible 
for this gene-specific recruitment of PRC2. To address this 
question, we performed ChIP experiments in the presence or 
absence of Trichostatin A (TSA) , a potent inhibitor of histone 
deacetylase activity. Mbd3 flox/ ~ ES cells were exposed to 
0.1 mM TSA for 2h. This treatment resulted in increased 
levels of H3K27 acetylation, which were comparable to 
those in Mbd3~ / ~ ES cells but which preceded gene expres- 
sion changes (Figure 5A; Supplementary Figure S5). In the 
presence of TSA, acetylation levels were increased at all 
genes tested, in contrast to the gene-specific increase in 
acetylation of H3K27 observed in Mbd3~ /_ ES cells. 
Exposure to TSA had little effect on the presence of Mi2p at 
previously identified target genes (Figure 5B) although Suzl2 
recruitment was reduced at all loci. This is unlikely to be 
a consequence of changes in protein abundance since levels 
of some NuRD component proteins decreased slightly in 
response to exposure to TSA, while PRC2 components 
showed a slight increase in abundance (Figure 5D). The 
observed dependency on acetylation levels for recruitment 
of PRC2 is consistent with previous reports that Suzl2 
occupancy shows an inverse correlation with the presence 
of H3K27ac at target genes (Pasini et al, 2010b). 

ES cells maintained for many passages in the absence of 
Mbd3 or PRC2 components may have undergone adaptive 
changes in response to prolonged culture. To determine 
whether PRC2 recruitment is immediately responsive to the 
presence of the NuRD complex, we made use of a tamoxifen- 
inducible Mbd3 system. The Mbd3b protein isoform was 



fused to the two copies of the mouse oestrogen receptor 
(MER-Mbd3b-MER), and expressed in Mb d3 -mill ES cells. 
In the absence of 4-hydroxytamoxifen, the protein is localised 
to the cytoplasm and the cells lack functional NuRD complex. 
Upon addition of tamoxifen, the protein translocates into the 
nucleus, restoring silencing of the NuRD/PRC2 target gene 
Html (Figure 6A and B). Recruitment of MER-Mbd3b-MER to 
the Html locus is observed by 20 h of tamoxifen treatment. 
In addition, over the same time period, H3K27 acetylation 
levels decrease and H3K27 trimethylation levels increase 
(Figure 6C). Coincident with these changes in H3K27 
modification and gene expression, there is a recruitment of 
Jarid2 to the Html locus. This experiment demonstrates that 
H3K27 modification status and PRC2 localisation at the Html 
locus are highly responsive to the presence or absence of a 
functional NuRD complex. Together with the lack of Suzl2 
binding at NuRD-specific genes in Mbd3~ /_ ES cells 
(Figure 4B), we conclude that NuRD-dependent deacetylation 
of H3K27 is required for PRC2 recruitment and subsequent 
trimethylation at NuRD target loci. 

Overlap of NuRD and PRC function in ES cells 

Both NuRD and PRC regulate gene expression and both are 
required for early embryonic development. However, the 
functions of these two complexes do not overlap completely 
as NuRD is required for lineage commitment of both early 
embryonic cells and ES cells, whereas ES cells or embryos 
lacking PRC2 components show some capacity for differen- 
tiation (Faust et al, 1995; O'Carroll et al, 2001; Pasini et al, 
2004; Shen et al, 2009; Leeb et al, 2010). In order to assess 
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Figure 5 Effect of acetylation on recruitment of PRC2 to chromatin. (A) ChIP for H3K27ac and IgG control in wild-type cells treated with TSA 
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with untreated. (C) Relative enrichment for Suzl2 and IgG control in wild-type cells treated with TSA compared with untreated. Error bars 
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Figure 6 Responsiveness of PRC2 recruitment to presence of 
NuRD. (A) Western blot of nuclear extracts prepared from wild- 
type, Aft>d3 _/ ~ and MER-Mbd3b-MER WSbd3~'~ cells after addition 
of 4-hydroxytamoxifen for indicated time. Alpha tubulin acts as a 
loading control. (B) Expression levels of Html in wild-type, 
Mbd3~ x ~ and MER-Mbd3b-MER Mbd3~'~ cells over the indicated 
4-hydroxytamoxifen time course relative to Mbd3 _/ ~ levels. (C) 
Occupancy levels for H3K27ac, H3K27me3, Jarid2 and Mer-Mbd3b- 
MER at the Html locus over the 4-hydroxytamoxifen time course in 
MER-Mbd3b-MER Mbd3~'~ ES cells. Levels are shown relative to 
highest point for each protein. 



more globally the degree to which NuRD- and PRC2- 
dependent transcriptional regulation overlaps in ES cells, we 
compared the gene expression changes found in Mbd3~ / ~ ES 
cells with those seen in ES cells lacking components of the PRCs 
(Leeb ef al, 2010). Using a significance cutoff of P<0.05, we are 
able to identify 1879 genes differentially expressed in Mbd3 _/ ~ 
ES cells also present on the microarray platform used in the 
Leeb et al study. Of these, ~ 23 % (436, Figure 7A, group I) were 
misregulated both in the Mbd3 and Polycomb mutant ES cells, 
consistent with the observed synergy in gene regulation be- 
tween the NuRD and Polycomb systems at a large proportion of 
target loci (Figure 7A). This subset of genes includes several 
encoding proteins important for embryonic development and 
key components of signalling pathways such as p-Catenin, 
Sfrpl, Tbx3, Tgfp and Wnt7b (Niwa et al, 2009; Pera and 
Tarn, 2010; Kelly et al, 2011; Supplementary Table 4). 

To identify more directly the extent to which PRC2 locali- 
sation is dependent on the action of NuRD, we performed 
ChlP-seq for Suzl2 in Mbd3~ / ~ ES cells. Of the 1928 high 
confidence Suzl2 binding peaks we identified in Mbd3 Flox/ ~ 
ES cells, 330 of these were not found in Mbd3~ f ~ ES cells 



(17.1%; Supplementary Table 5) . Consistent with being PRC2 
targets, these genes were generally enriched for trimethyla- 
tion and depleted for acetylation at H3K27 (Supplementary 
Figure S6) . In contrast, this same set of genes show a relative 
decrease in H3K27Me3 levels and a corresponding increase in 
H3K27Ac levels in Mbd3-null ES cells. Notably, the loss of 
methylation is most pronounced close to the transcription 
start sites (Supplementary Figure S6). These data strongly 
support a model whereby NuRD and PRC2 control expression 
through a balance between acetylation and methylation 
status of H3K27. 

Of the genes displaying NuRD-dependent peaks of Suzl2 
binding, 133 were also identified as Mi2p targets by ChlP-seq 
(Table I). This group is highly enriched for genes, which have 
functions critical for transcriptional regulation (Table II). 
Many of these transcription factors have known functions 
in stem cell maintenance or pluripotency (e.g., Tbx3, Klf4 
and Foxd3; Hanna et al, 2002; Niwa et al, 2009) or early 
embryogenesis and differentiation processes (e.g., Gata5, 
Sfrpl and Smarcd3; Lessard et al, 2007; Lou et al, 2011). 
Thus, while the number and identity of genes for which 
NuRD function appears to be required for PRC2 binding 
represents a relatively small proportion of either total NuRD 
or PRC2 targets, many of the genes targeted for joint regula- 
tion by both PRC2 and NuRD are known to be critical 
for development and highlights the importance of such a 
mechanism in vivo. Based upon these findings, we conclude 
that the mechanism described here underlies the regulation 
of a significant subset of PRC and NuRD target genes in ES 
cells, and that this interaction may be responsible for the 
shared biological functions of the two silencing complexes. 

Discussion 

NuRD-mediated gene regulation has been shown to facilitate 
developmental transitions in early mouse embryos, lineage 
commitment in ES cells and developmental decisions in 
haematopoietic and epithelial stem cells (Kaji et al, 2006, 
2007; Kashiwagi et al, 2007; Yoshida et al, 2008). Using a 
combination of expression analysis and ChIP experiments, 
we demonstrate for the first time a mechanism by which 
NuRD regulates the transcription of specific genes in ES cells. 
We show that NuRD-dependent deacetylation of H3K27 
makes NuRD target genes available for further repressive 
action by the PRC2. Our data support a model for a two- 
step process for transcriptional silencing of bivalent genes 
(Figure 7B) centering on a switch between the acetylation 
and methylation status of H3K27. 

Our model is consistent with previous reports that describe 
both reciprocity between H3K27 acetylation and methylation 
and correlation of H3K27 modification status with transcrip- 
tion (Tie et al, 2009; Jung et al, 2010; Pasini et al, 2010b). The 
status of H3K27 can be reversibly modified by a combination 
of enzyme activities (Figure 7B; Tie et al, 2009; Pasini et al, 
2010b), allowing an exquisitely responsive yet robust me- 
chanism of gene regulation. Repressive modification (di- and 
trimethylation) is carried out at H3K27 by histone methyl- 
transferase activity (i.e., by Ezh2 in the PRC2 complex) and 
removed by the action of histone demethylases (such as 
tri thorax group proteins). Histone H3K27 may be acetylated 
to establish active transcription (by histone acetyltrans- 
ferases) or deacetylated (by histone deacetylases) as a further 
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level of regulation. The identity of some of these enzymes has 
been determined (Tie et al, 2009; Pasini et al, 2010b). 
Regulation of transcriptional status in the case of H3K27 



Table I Mi2|3 target genes showing Mbd3 -dependent Suzl2 binding 





IV 



B 



Deacetylation 
(gene specific) 



Demethylation 




NuRD 




H3K27ac 



H3K27 



H3K27me3 





PRC2 



Acetylation 



TRANSCRIPTION 



Methylation 
(all bivalent genes) 




TRANSCRIPTION 



Demethylation 




H3K27ac 



H3K27 H3K27me3 
PRC2 




TRANSCRIPTION 



1700029F12Rik 


E030030I06Rik 


Klf4 


Pthlr 


1700067P10Rik 


Edn2 


Lhx5 


Pthlh 


2810030E01Rik 


Efna2 


Lingo 1 


Ptprd 


6030419C18Rik 


Egr4 


Lphn2 


Ptrf 


8430427H17Rik 


Ezr 


Ly6c2 


Raplgap2 


9030425EllRik 


Faml29c 


Lyzl4 


Reml 


A930004D18Rik 


Fam92b 


Manlcl 


Rnf220 


Abcc5 


Fosl2 


Mapkl5 


Ror2 


Abhdl5 


Foxd3 


Mapt 


RP23-197C15.2 


AC 131 780. 10 


Gata5 


Mcf21 


Rufy4 


AC131 780.6 


Gbxl 


Mlxipl 


Scarf2 


AC140264.1 


Gbx2 


Nav2 


Sdk2 


AC140331.1 


Gm5089 


Ncaml 


Sfrpl 


Agpat9 


Gm5607 


Ncoa2 


Slc22a21 


AL837506.1 


Gm5627 


Nebl 


Smarcd3 


Ankrd56 


Gm6020 


Nfil3 


Sox21 


Ano4 


Gm6304 


Nkx6-2 


Tbca 


Astn2 


Gpm6b 


Nmnat2 


Tbx3 


Atxnl 


Gprl50 


Npas2 


Tcf4 


Bach2 


Gpr83 


Nr2f2 


Tcte2 


Bahccl 


Gramdlb 


Ntnl 


Thpo 


Brunol4 


GrblO 


01fm2 


Tmtcl 


Camtal 


Grik3 


01ig3 


Unc5b 


Caszl 


Grik4 


Ovoll 


Unc5d 


Cbfa2t3 


Grin2a 


P4ha2 


ZbtblO 


Cbln4 


Hmga2 


Pafld 


Zbtb7c 


LQko 


Hrh2 


Pcangall 


Zeb2 


Clstnl 


Hs3st3bl 


Pde8a 


Zfhx3 


Cnih3 


Htral 


Pdgfa 


Zfp3612 


Crem 


Igsf21 


Pou2fl 


Zfp423 


Crlfl 


Irx2 


Ppp2r2c 


Zfp804a 


Csf2ra 


Kazaldl 


Prickle 1 




Cuxl 


Kcnh3 


Prkarlb 




Dgkz 


Kcnk3 


Prkcdbp 





Figure 7 Overlap of NuRD and PRC function in gene regulation. 

(A) Expression patterns of genes differentially expressed between 
wild-type and Mbd3~ / ~ ES cells (P<0.05) were assessed in data 
from Leeb et al (2010) describing gene expression changes in PRC 
mutant ES cells. A total of 2907 probe sets were found to be 
differentially expressed between the reference and Mbd3 mutant 
samples, mapping to 1967 unique genes using the re-annotated 
array information provided by Leeb et al. Of these, 1879 genes were 
also present on the Leeb microarray platform, although not neces- 
sarily differentially expressed. The log 2 fold change of these genes 
relative to their respective wild-type samples is illustrated in the 
figure. Genes showing reduced expression compared with wild-type 
cells are indicated in blue; those showing increased expression are 
indicated in yellow. Mbd3_l, Mbd3_2 and Mbd3_3 indicate results 
from three different Mbd3~ / ~ ES cell samples corresponding to two 
independently derived ES cell lines. RinglB_l-3, Eed_l-3 and 
RinglB/Eed_l-3 represent results from three replicates each of 
Ring IB-null ES cells, Eed-null ES cells and RinglB/Eed-double 
null ES cells, respectively (Leeb et al, 2010). Where multiple probes 
map to the same gene, only the probe displaying the largest 
deviation from the reference sample is included here. The four 
main gene clusters are indicated on the right hand side (I-IV). 

(B) Two-step model for control of transcription illustrating relation- 
ship between acetylation status, here controlled by NuRD, and 
methylation of H3K27 by PRC2. Top panel: in wild-type cells, 
PRC2 and NuRD function fully. NuRD directs deacetylation of 
specific genes, which then become available for trimethylation by 
PRC2. Middle panel: in cells lacking functional PRC2, for example, 
EedT^ cells, H3K27 trimethylation is lost genome-wide, the balance 
moves towards the acetylation/deacetylation cycle and transcrip- 
tion of bivalent genes increases overall. Bottom panel: in the 
absence of NuRD, deacetylation fails to occur at specific loci only. 
At these genes, there is an increase in acetylation, a subsequent 
reduction in trimethylation at H3K27 through loss of substrate for 
PRC2 and an increase in transcription. 



©2012 European Molecular Biology Organization 



The EM BO Journal VOL 31 | NO 3 | 2012 601 



NuRD and PRC2 interact to control gene expression 

N Reynolds et al 



Table II Gene ontology analysis of NuRD-dependent Suzl2 targets 



GOMFID 


P-value 


Count 3 


Size b 


Term 


GO:0003700 


<0.01 


23 


738 


Sequence-specific DNA binding transcription factor activity 


GO:0043565 


<0.01 


14 


430 


Sequence-specific DNA binding 


GO:0030528 


<0.01 


11 


286 


Transcription regulator activity 


GO:0016564 


<0.01 


9 


206 


Transcription repressor activity 


GO:0000975 


<0.01 


9 


209 


Regulatory region DNA binding 


GO:0016563 


<0.01 


9 


335 


Transcription activator activity 


GO:0010843 


<0.01 


7 


189 


Promoter binding 


GO:0008270 


<0.01 


18 


1358 


Zinc ion binding 


GO:0005488 


<0.01 


18 


2712 


Binding 


GO:0005515 


<0.01 


23 


3229 


Protein binding 


GO:0005216 


<0.01 


7 


367 


Ion channel activity 


GO:0015267 


<0.01 


7 


387 


Channel activity 


GO:0003677 


<0.01 


11 


1076 


DNA binding 


GO:0060089 


<0.01 


23 


2850 


Molecular transducer activity 


GO:0003690 


<0.01 


4 


140 


Double-stranded DNA binding 


GO:0003702 


<0.01 


4 


148 


RNA polymerase II transcription factor activity 


GO:0005102 


<0.01 


10 


877 


Receptor binding 


GO:0003705 


0.01 


3 


82 


Sequence-specific enhancer binding RNA polymerase II transcription factor activity 


GO:0035257 


0.01 


3 


82 


Nuclear hormone receptor binding 


GO:0008022 


0.01 


4 


169 


Protein C-terminus binding 


GO:0030165 


0.01 


3 


94 


PDZ domain binding 


GO:0016566 


0.01 


3 


97 


Specific transcriptional repressor activity 


GO:0008134 


0.01 


5 


299 


Transcription factor binding 


GO:0008201 


0.01 


3 


104 


Heparin binding 


GO:0046872 


0.01 


23 


3246 


Metal ion binding 


GO:0043167 


0.01 


23 


3297 


Ion binding 


GO:000S179 


0.01 


3 


118 


Hormone activity 


GO:0003682 


0.01 


4 


217 


Chromatin binding 



a Number of genes from the test set [n = 133) corresponding to this GO term. Only those categories containing three or more genes, and with a 
P-value <0.01 are shown. 

b Number of total genes corresponding to this GO term. 



may therefore be regarded as two interconnected cycles, with 
methylation/demethylation establishing the silencing of a 
gene, and acetylation/deacetylation controlling full activation 
of transcription (Figure 7B) . 

To date, models for transcriptional regulation revolving 
around H3K27 status have lacked any mechanism for locus 
specificity. For instance, loss of Suzl2 in ES cells results in 
loss of H3K27 trimethylation, and an increase in H3K27 
acetylation on a global scale (Pasini et al, 2010b). Similarly, 
knockdown experiments to deplete cells of the histone acet- 
yltransferases Cbp (Crebbp) or p300 (Ep300) resulted in 
changes in H3K27 acetylation levels that were detectable in 
western blots of bulk histones (Pasini et al, 2010b) . This view 
is in contrast to the activity of NuRD described here, in which 
transcription and H3K27 acetylation of a relatively small 
number of genes are affected (Figure 2). Here, we show 
that the balance between acetylation and methylation status 
of H3K27 can be controlled at a fine level by NuRD, which 
indirectly affects the activity of PRC2 by controlling substrate 
availability at NuRD target genes. This presumably takes the 
form of making unmodified H3K27 available for monomethy- 
lation prior to di- and trimethylation by the PRC2 complex. 

Based on our model (Figure 7B), several mechanistic 
predictions are possible. In the event of PRC2 being removed 
from the system (Figure 7B, middle panel), H3K27me3 will be 
lost and H3K27 acetylation will increase globally, a fact that 
has been established in other studies (Tie et al, 2009; Jung 
et al, 2010; Pasini et al, 2010b) and confirmed here 
(Figure 2B). The resulting shift from a silent methyl mark 
to that of acetylation is accompanied by a general upregula- 
tion of bivalent genes in the absence of PRC2. A lack of 
methylation at H3K27 would not necessarily affect recruit- 
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ment of NuRD to the chromatin as the acetylation/deacetyla- 
tion element of the cycle remains in balance and substrate 
levels will be unaffected (Figure 7B), an hypothesis that we 
confirmed experimentally for the NuRD component Mi2p 
(Figure 4C, E). 

As NuRD acts in a gene-specific manner, its loss affects 
acetylation levels of H3K27 at direct target genes only. Here, 
accumulation of H3K27 in an acetylated state would result in 
localised loss of substrate for PRC2 and subsequent reduction 
in occupancy only at specific gene loci (Figures 4B, D and 7B, 
bottom panel). We propose that deacetylation of H3K27 is 
required for PRC2 function, and that the action of NuRD 
enables PRC2 to silence specific loci. Pasini et al (2010b) 
proposed that the function of PcG proteins is to prevent 
acetylation of H3K27. According to our model, the action of 
enzymes which affect the methylation status of this histone 
residue will inevitably affect the availability of substrate for 
enzymes which control the acetylation status, and vice versa. 
Since it is the balance between the activation cycle (acetyla- 
tion/deacetylation) and the silencing cycle (methylation/ 
demethylation) that ultimately determine transcriptional 
status, an increase in influence of either one will necessarily 
have an affect on the other. 

While a requirement for NuRD in ES cell differentiation 
has been established (Kaji et al, 2006), the identities of 
specific gene loci targeted by the complex have remained 
largely unexplored. In this study, we have shown a correla- 
tion between promoter occupancy and gene expression levels 
using a combination of ChIP and microarray analysis in 
Mbd3 _/_ ES cells. Based on our expression results, there 
appears to be no bias towards target genes with H3K4me3, 
H3K27me3 or bivalent histone modifications. In fact, while 
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NuRD is traditionally viewed as a transcriptional co-repres- 
sor, the majority of genes expressed in its absence exhibit 
regions of H3K4me3 at their promoter regions, indicative of 
actively transcribed genes. This is consistent with a recent 
study in which a number of different histone deacetylases 
were found to associate globally with actively transcribed 
genes in human T cells (Wang et al, 2009). Therefore, rather 
than acting as a conventional silencer, our data are consistent 
with NuRD functioning as a modulator of transcriptional 
activity at a large proportion of its target genes. Indeed, we 
have recently shown that NuRD controls the dynamic range of 
expression levels of some actively transcribed genes in ES cells 
to facilitate lineage commitment (Reynolds et al, submitted) . 

Interactions between the Mi2 and PcG family proteins have 
been reported to act either in parallel with each other, for 
example, at sites of DNA damage (Chou et al, 2010), antag- 
onistically, for example, in the control of cell identity in plants 
(Aichinger et al, 2009), or in synergy with each other, as in 
the case of Drosophila dMi-2 which has been shown to 
specify sites of Polycomb-mediated repression in vivo 
(Kehle et al, 1998). How these interactions result in transcrip- 
tional changes has not been determined. Using gene expres- 
sion profiling in ES cells, we show the extent of functional 
overlap of the NuRD- and PRC-dependent gene expression 
programmes. We further provide the first clear evidence for a 
mechanism by which these important regulators may interact 
to control gene expression in mammalian cells, explaining 
observations made previously in both plants and human 
leukaemic cells (Morey et al, 2008; Zhang et al, 2008a). 
Uniquely, we show that NuRD targets the action of the 
PRC2 complex to specific gene loci via H3H27 deacetylase 
activity. We further show that NuRD activity is required 
for the recruitment of PRC2 to a set of developmentally 
important transcriptional regulators in ES cells, providing a 
mechanism for targeting PRC2 activity to specific loci. In 
summary, here we provide evidence that NuRD and PcG act 
together to both dictate and reinforce transcriptional 
silencing in ES cells. 

Materials and methods 

ES cell culture 

ES cells were grown in self-renewing conditions, that is, in the 
presence of serum and murine LIF. Mbd3-null and Eed-null ES cell 
lines have been described previously (Silva et al, 2003; Kaji et al, 
2006). Eed-null ES cells were cultured using a layer of feeder cells 
which were depleted from ES cells upon harvest by taking 
advantage of their differential attachment properties and allowing 
the feeders to attach to non-gelatinised plates. TSA (Sigma) was 
added to normal ES cell media to a final concentration of 0.1 mM for 
2 h prior to harvesting and further analysis. An expression construct 
in which the Mbd3b coding region was fused at both N- and 
C-termini to mouse oestrogen receptor (MER) domains (Verrou 
et al, 1999) was transfected into ES cells. Clones in which 

uninduced cells showed neither nuclear localisation of the MER- 
Mbd3b-MER protein by immunofluorescence nor silencing of 
known NuRD target genes were selected for further analysis. 
Responsiveness to tamoxifen was measured both by immunofluor- 
escence and quantitative RT-PCR for NuRD responsive genes. 

RNA isolation and gene expression analysis 

Total RNA was prepared from Mfad3-null or parental ES cell lines 
using Trizol (Life Technologies) according to the manufacturer's 
protocol and hybridised to Mouse WG-6 Expression arrays 
(Illumina, Inc). Fluorescence data were processed using software 
packages from v2.7 of the Bioconductor project (Gentleman et al, 
2004). Data files were handled using beadarray v2.0.0 (Dunning 



et al, 2007) and hybridisation quality assessed with arrayQuality- 
Metrics v3.2.0 (Kauffmann et al, 2009). Normalisation was 
performed using the ssn algorithm from lumi v2.2.0 (Du et al, 
2008). Microarray probe sets with an interquartile range <0.4 were 
removed, and for those remaining differential expression was 
calculated using limma v3.6.0 (Smyth, 2005). The resulting 
P-values were corrected for multiple testing using the Benjamini 
and Hochberg (1995) approach, and probes with P<0.01 were 
deemed significant. To improve the interpretation of the data, 
microarray probe identifiers were re-annotated as outlined in 
Barbosa-Morais et al (2010). Published primary microarray data 
from EedT^ ES cells (Leeb et al, 2010) were re-analysed using RMA 
(Irizarry et al, 2003) prior to comparison with data obtained in this 
study from Mbd3^ ES cells. 

Quantitative real-time PCR 

Gene expression levels were verified by quantitative PCR carried 
out in triplicate using cDNA from total RNA prepared from ES cell 
lines with Fast SYBR green Master Mix (Life Technologies) on the 
7900HT Fast Real-Time PCR System (Life Technologies). Primers to 
either fi-actin or PplA were used as endogenous controls and 
relative values calculated using the ACt method. Primer sequences 
are listed in Supplementary Table 6. 

Chromatin immunoprecipitation 

ChIP was carried out using antibodies to endogenous proteins listed 
in Supplementary Table 7 and standard ChIP protocols. In brief, 
cells were fixed using 1% formaldehyde for lOmin at room 
temperature and fixation quenched with 150mM glycine. For 
Mi2[i ChlP-seq, formaldehyde fixation was preceded by 45min 
incubation in 2mM DSG (disuccinimidyl glutarate (Sigma); Wang 
et al, 2009). Chromatin was sheared by sonication to an average 
fragment size between 200 and 300 bp using a Bioruptor sonication 
instrument (Diagenode). Immunoprecipitation was performed 
using Protein A sepharose, pre-blocked with fish gelatin and 
single-stranded herring sperm DNA. Quantitative PCR of ChIP 
DNA was carried out as above using gene-specific primers listed in 
Supplementary Table 6. All qPCR was carried out in triplicate for at 
least three biological replicates. Analysis was carried out either as 
percentage of input DNA or as ACt values in order to compare 
relative levels between samples, using primers to the promoters of 
Actin or Cdx2 (Supplementary Table 6b) as endogenous controls. 
P-values were calculated by T-test analysis. 

ChlP-seq analysis 

For sequencing of ChIP DNA, samples from three individual ChIP 
experiments for each cell line were verified by qPCR before pooling for 
library construction. Sequencing was carried out on the Illumina GAIIx 
at a read length of 36 bp. These were mapped to build mm9 of the 
mouse reference genome using vO.12.7 of the Bowtie short-read aligner 
(Langmead, 2010) . Redundant reads were removed, and those from the 
two biological replicates were merged into a single data set. MACS 
(Zhang et al, 2008b) and CCAT (Xu et al, 2010) peak-finding programs 
were used to identify histone-modified regions by comparing both 
Mbd3-null and parental samples to their respective input DNA controls. 
MACS was used to identify peaks in Mi2(i and Suzl2 ChlP-seq data. 
Equal numbers of reads from ChIP and control samples were used for 
each comparison, and regions that were determined to be significantly 
enriched by both peak callers were included in downstream analyses. 
Annotation of modified regions was performed with PeakAnalyzer 
(Salmon-Divon et al, 2010) and CEAS (Shin et al, 2009). 

The Refseq gene model annotation was used to define transcrip- 
tion start and end sites. The CEAS tool (Shin et al, 2009) was used 
to extract the signal surrounding every relevant genomic region. 
A region comprising three kilobases upstream and downstream of 
each feature was divided into 120 bins (50 bp each). The median 
signal of each bin was calculated, and then averaged across all 
features. The resulting average profiles were normalised to input 
signal level. ChlP-seq and expression data have been deposited in 
the ArrayExpress database under accession numbers E-MTAB-888 
and E-MTAB-889, respectively. 

Western blotting 

Nuclear extracts were prepared as described previously (Kaji et al, 
2006), subjected to SDS-PAGE and probed with antibodies listed in 
Supplementary Table 4. 

Acid extraction of bulk histones was achieved by washing nuclei 
in 0.1 M HC1 according to the standard protocols. 
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Supplementary data 

Supplementary data are available at The EMBO Journal Online 
(http://www.embojournal.org) . 
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